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Absolute  instability  of  G  rtler  vortices  on  the  severely  curved  concave  pressure 
side  of  a  gas-turbine  blade  is  the  main  thrust  of  the  third-year  work  under  this  grant.  For 
the  most  part,  the  G  rtler  vortices  have  been  investigated  in  an  incompressible  boundary 
layer  over  thin-wing  sections  or  artificial  inserts  on  an  otherwise  flat  plate.  The  cascade 
of  modem  aircraft  engines  operate  in  the  high  subsonic  Mach  number  regime  with 
velocity  fields  strongly  affected  by  centrifugal  forces  maintained  by  the  large  curvature  of 
profiles.  Unsteady  spiral-type  vortices  developing  in  these  environments  provoke  the 
absolute  instability  in  the  streamwise  direction  of  the  boundary  layer  leading  to  earlier 
transition.  An  effort  undertaken  after  the  meeting  in  Shalimar  ( May  29-31,  2002  )  show 
that  the  heat  transfer  coefficient  is  even  more  susceptible  to  enhancing  oscillations  in  the 
upstream  moving  wave  packets  than  the  pressure. 

The  essence  of  the  matter 


The  shape  of  a  blade  in  broad  use  in  the  modem  gas-turbine  engines  is  shown  in 
Figure  1.  The  concave  pressure  side  is  exposed  to  severe  loading  that  can  cause  bad 
damages  (spallation  of  thermal  barrier  coatings).  The  large  curvature  of  the  concave 
pressure  side  maintains  centrifugal  forces  giving  rise  to  G  rtler  vortices  which  make  the 
problem  all  the  more  acute.  The  interaction  of  Tollmien-Schlichting  waves  and  unsteady 
spiral-type  vortical  disturbances  results  in  upstream-advancing  wave  packets  that  are  at 
the  heart  of  absolute  instability  in  the  streamwise  direction  of  the  oncoming  flow.  This 
kind  of  instability  was  discovered  by  Ryzhov  (2001)  in  incompressible  boundary  layers 
and  recently  proved  by  Ryzhov  &  Bogdanova-Ryzhova  (2002)  to  exist  in  the  high 
subsonic  Mach  number  regime.  Surprisingly  enough,  temperature  pulsations  in  the 
upstream-moving  wave  packets  appear  to  have  even  larger  amplitude  than  the  pressure  if 
the  characteristic  length  of  G  rtler  vortices  grows  without  bound.  The  heat  transfer 
coefficient  also  increases  with  the  vortex  length  (  Ryzhov  &  Bogdanova-Ryzhova,  2003). 
The  temperature  disturbances  are  triggered  by  the  pressure  and  velocity  fields  that  need 
to  be  evaluated  first. 

Conventional  View  of  the  G  rtler  instability  mechanism 

Typically,  the  G  rtler  vortices  are  thought  of  as  steady  disturbances  having  no 
characteristic  length  in  the  streamwise  direction  (  Saric,  1994).  Figure  2  illustrates  this 
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conventional  point  of  view  under  the  tacit  assumption  of  incompressible  flow. 
Theoretically,  an  appropriate  representation  of  velocity  components  reads 

t{x,y,z)  =  L{x,yj  +  r/(y,)eflxCOs(27tz/Xz)  (la) 

v(- x,y,z) )  =  V(x,y)  +  ^(yje^  COs[lnz/X  z)  (lb) 

v\x,y,z)  =  w(j)e^sin(27iz/^z)  (lc) 

Here,  the  real  exponent  P  determines  the  disturbance  amplification  rate  in  the 
streamwise  direction,  but  there  is  no  reference  wavelength  in  this  direction.  The  reference 
wavelength,  X  z ,  is  specified  only  in  the  spanwise  direction.  It  can  be  shown  that  stability 
properties  of  the  boundary  layer  on  a  concave  wall  depend  on  a  parameter 

G  =  k  Re  (2) 

called  the  G  rtler  number  that  incorporates  the  surface  curvature  k  and  the  Reynolds 
number  Re>>  1.  A  number  of  other  parameters  enter  the  problem  posed  for 
compressible  flow.  The  steady  G  rtler  vortices  do  not  disturb  the  boundary  layer 
upstream  of  an  external  agency. 

Unsteady  G  rtler  vortices 

Unsteady  spiral-type  vortical  disturbances  have  been  recorded  many  years  ago  in 
a  related  ( incompressible  )  flow  in  a  gap  between  two  co-axial  cylinders,  one  of  which 
rotates  and  the  other  is  at  rest.  The  spiral-type  vortices  create  a  wavy  motion  shown  in 
Figure  3  as  applied  to  the  problem  of  interest  on  the  boundary  layer  developing  on  a 
concave  solid  surface.  Here  X  stands  for  the  wavelength  in  the  spanwise  direction, 
however  a  characteristic  wavelength  is  clearly  seen  to  form  in  the  disturbance  pattern  also 
in  the  direction  of  the  oncoming  stream. 

The  problem  has  been  attacked  theoretically  by  Denier,  Hall  &  Seddougui  (1991) 
who  identified  five  different  regimes  inherent  in  the  high-Reynolds  number  flow.  Starting 
from  the  results  of  this  paper,  Choudhari,  Hall  &  Streett  (1994),  arrived  at  a  wrong 
conclusion  that  no  upstream  influence  results  from  the  viscous  /  inviscid  interaction 
controlled  by  centrifugal  forces.  As  a  matter  of  fact,  the  first  experimental  evidence  of 
upstream  influence  exerted  by  the  vortical  disturbances  on  a  concave  surface  can  be 
traced  back  to  much  earlier  work  by  Mangalam,  Dagenhart,  Hepner  &  Meyers  (1985). 

An  extended  version  of  the  triple-deck  scheme 

An  asymptotic  approach  based  on  the  assumption  that  the  Reynolds  number  takes 
on  sufficiently  large  values  is  required  to  resolve  this  controversy.  The  triple-deck  theory 
offers  a  clue  to  start  mathematical  treatment.  However,  the  classical  version  of  the  theory 
is  based  on  the  scaling  intrinsic  to  Tollmien-Schlichting  waves  (  Smith,  1979;  Zhuk  & 
Ryzhov,  1980  ).  Vortical  eigenmodes  are  not  covered  by  this  scaling  and  therefore  fall 
beyond  the  Tollmien-Schlichting  spectrum  of  frequencies.  Neglect  of  centrifugal  forces 
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maintained  by  the  surface  curvature  is  the  decisive  factor  missing  from  the  classical 
triple-deck  analysis  of  wave  motion.  Higher  order  terms  have  to  be  retained  in  the 
asymptotic  expansions  to  account  for  the  curvature  effects  giving  rise  to  G  rtler  vortical 
eigenmodes. 

Guided  by  this  preliminary  consideration  we  introduce  the  time-space  scaling 

f=  z2  ~j~t,  **=  L*(l+  z2x),  2*=  e 2L*z  (3) 

^00 

typical  of  the  Tollmien-Schlichting  spectral  range  fixed  in  terms  of  a  small  parameter 

l 

e  =  Re  *  (4) 

for  Re>>  1.  The  definition  of  the  normal-to-wall  distance  depends  on  a  specific 
sublayer. 

The  Prandtl  variable  y2  furnishes  a  pertinent  scaling 


y*=  zAL*y2 

(5) 

for  most  of  the  boundary  layer  where  the  desired  functions  have  the  following  asymptotic 
expansions 

*  =  ^ ^21  E  ^24  “ " " 

^00 

(6a) 

*  =  V2  =  8  ^21  (^^^2  amm^  ^  ^24  3^2 5  "  ’ 

Mdo 

(6b) 

~“t~  =  w2  =  z2w2l(t9x,y2,z)+*.Az5w24(t,x,y2,z)+... 

^00 

(6c) 

P *  K2  -  Pi  =  z2p2](t,x,y2,z)+...+  zAp23(t,x,y2,z)  +  z5p24(t,x,y2,z)+... 

P  00^00 

(6d) 

lir 

~  P2  =  ^P2l(^5**'53;2  P  24(^5  ^y2  ■ 

P  00 

(6e) 

The  initial  velocity  and  density  profiles  in  the  unperturbed  boundary  layer  on  the 
concave  side  of  a  blade  are  given  by  the  leading-order  terms  U0[y2)  and  ^(j^)  in  (6a) 
and  (6e),  respectively.  Logarithmic  terms  are  not  explicitly  indicated  since  they  have 
nothing  to  do  with  the  instability  problem  in  question.  The  scaling  of  the  velocity, 
pressure  and  density  fields  comes  from  the  classical  triple-deck  analysis  providing  the 
opening  stage  of  the  present  study. 

Based  on  the  ideas  set  forth  at  the  beginning  of  this  report,  let  us  focus  on  the 
pressure  variations  across  the  boundary  layer  which  are  balanced  out  by  centrifugal 
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forces  supported  by  the  blade  surface  bending.  Let  designate  the  contribution  to  the 
pressure  from  the  curvature  effects.  Then 


Ak)  =  e4 


x,y2)-ni^t(t,x,y2,z) 


+  .  .  . 


(7) 


According  to  the  classical  boundary-layer  theory  by  Prandtl, 

?MK)  «  2 

Sy2  ‘  '  ViteP^ 

=  -s  K  | (jPh )  ^  (.^2 )p 2 1  (^> ^>y2 > ^)]"^  •  ■  * j 


Hence 

=  -K^)U)^U)  (9) 

is  the  pressure  gradient  in  the  initially  unperturbed  flow  and 
a  M 

~~~=  -«[2^Wl6(>4Hi(/‘»-*04»z)  +  d°) 

determines  the  self-induced  pressure  oscillations  responsible  for  the  boundary-layer 
instability  provoked  by  centrifugal  forces. 

To  evaluate  the  last  term  giving  the  d{z 5 )  -contribution  to  the  pressure  we  need 
only  the  first-order  solution  for  velocity  components  and  density.  This  solution  can  be 


explicitly  expressed  in  the  form 

%h\  =  ^ 

(11a) 

(lib) 

2'  ■  -  f; 

(11c) 

a  \dR§ 
p21  -  At,x,z) 

(lid) 

where  -  Aj,  x,  z)  has  a  simple  meaning  of  the  instantaneous  displacement  thickness. 
Using  (1  la,d),  the  self-induced  pressure  becomes 
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(12) 


dy2  dy2 


Integration  of  (12)  allows  the  thermal  conditions  on  the  blade  surface  to  be  taken  into 
account  when  evaluating  the  heat-transfer-coefficient  disturbances.  The  final  result 
implicitly  depends  on  the  Mach  number  involved  in  the  velocity  and  density 

distributions,  LJ^y2 )  and  ^(y2) ,  across  the  boundary  layer. 


Interaction  law 

As  (12)  shows,  the  self-induced  pressure  oscillations  directly  relate  to  the 
instantaneous  displacement  thickness  -  Ay,  x,  z) .  The  interaction  law  to  connect  these 
two  quantities  derives  from  the  study  of  the  disturbance  pattern  in  the  outer,  essentially 
inviscid  sublayer  where  the  normal-to-wall  distance  is  scaled  as 

y*=  e3L*yx 

We  start  again  with  the  asymptotic  expansions 
— ~r=  Wj  =  1+  E2t*lx(t,x,yx,z)+...+  E5uX4(t,x,yl,z)+... 

~TF~  Vl  = 

Mdo 

—r  =  w\  -  z2wxx(t,x,yx,z)+  ...Jr£5wX4(t,x,yx,z)Jr ... 

★  * 


*  ±2  ”  Pi  “  £  •  ••+  8  p^4(t9x9y^9zj+ 

PM 

(14d) 

=  p  \  =  1  +  £  p11(^^lJz)+...+  £  p  .«■ 

(14e) 

characteristic  of  the  Tollmien-Schlichting  spectral  range.  In  normalized  variables 

{x,z}  =  (l-  Adpx 
_  1  1  1 

P=(  1-A^)  **IPP 

A=(l- 

the  first-order  solution 


(15a) 

(15b) 

(15c) 


(13) 

(14a) 

(14b) 

(14c) 
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a2^<,u)M2 


(16) 


1  1  00  CO 

P\X(t,x,0,z)=  (l-  At)  2  f 


(l -h£)  \x-^)2  +  (z- qf 


incorporating  an  explicit  dependence  on  the  Mach  number  equally  applies  to  flat  as  well 
as  curved  surfaces.  This  should  be  supplemented  with  a  term  t3 proportional 


to  the  surface  curvature  k  .  Making  use  of  the  results  of  the  preceding  section  yields  an 
extended  version  (Ryzhov,  2001;  Ryzhov  &  Bogdanova-Ryzhova,  2002) 


1  1  OO  00 

£=--(i-A£pteJ  — 

*£)  '(*- 02  +  (4_  ?)2 

D=  (l  -A£fTjd{Z/£y 


—dq  +  s  k  DA 


(17a) 


(17b) 


of  the  interaction  law  where  and 

_3 

S=  k  Re  8  D  (18) 

are  two  similarity  parameters.  It  is  worth  noting  that  the  generalized  interaction  law 
consists  of  two  components.  The  first  term  on  the  right-hand  side  of  (17a)  is  associated 
with  the  Tollmien-Schlichting  waves,  the  second  one  relates  to  the  spiral-type  G  rtler 
vortices  induced  by  centrifugal  forces.  Thus  (17a,b)  stem  from  the  interaction  of  the 
wave  and  vortex  eigenmodes  evolving  on  the  concave  pressure  side  of  a  turbine  blade. 

Final  formulation 


The  interaction  law  (17a,b)  defines  both  components  dp/dx  and  dp/oz  of  the 
pressure  gradient  in  the  Prandtl  equations 


du  dv  dw 
—  +  — +  —  =  0 
d%  dp  82 

8ti  du  du 


dw  dw  ^  dw  ^  dw 

~di  +u~d&  +  vH$>  +  m1)2 


dp  d2u 
djc+  dp2 
dp  d2w 
~  dz+  dtf 


(19a) 

(19b) 

(19c) 


ensuing  from  the  system  of  original  Navier-Stokes  equations  in  the  viscous  near-wall 
sublayer.  This  is  a  crucial  simplification  of  the  problem  in  the  framework  of  the 
asymptotic  high-Reynolds  number  formulation  based  on  rescaled  normalized  variables 
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1  3  1 


L'/l/. 

y 

V 


Uw 


'rr\ 

(20a) 

~  $ 

13  5  3 

85(l-A£pt  JO(Tjiyy 

(20b) 

Jj\  =  e(l-  a£)  *xAwd{fJ£)2{u,v)i 

(20c) 

13  3  1 

e3(l  -AtYxld{TjtYv 

(20d) 

the  notation  being  standard.  The  interaction  law  preserves  the  ellipticity  of  the  original 
Navier-Stokes  equations  which  brings  about  upstream  influence.  Mathematically,  it  is 
just  this  property  which  underlies  absolute  instability  in  the  streamwise  direction.  The 
limit  conditions 


v*  \  dw  1  dp 
u- y->  At,x,z),  as  $>->  oo 

complete  the  asymptotic  formulation  for  Re  >>  1 . 

A  comparison  of  classical  and  interactive  approaches 


(21a,b) 


The  G  rtler  number  defined  in  (2)  gives  the  similarity  parameter  for 
incompressible  flow  in  the  scope  of  classical  hydrodynamic  stability  theory.  The  skin 
friction  xw varies  along  the  vortex  length.  The  interactive  boundary  layer  scheme  predicts 
a  different  similarity  parameter 


k  Re~*(l-  A£)*xjd{fJ£f 


(22) 


3 

depending  on  k  Re  8  rather  than  k  Re  and  involving  a  combination  of  the  Mach  number 
characteristic  skin  friction  xw.  Chapman  viscosity  constant  C  and  temperature 

ratio  Twj Tw .  All  these  additional  parameters  should  be  considered  separately  in  the 

classical  approach  rendering  computed  results  hardly  conceivable.  Thus,  both  the  scaling 
and  similarity  parameters  prove  to  be  quite  different.  Hence,  the  unsteady  spiral-type 
vortices  are  not  an  extension  of  the  classical  G  rtler  vortex  structures,  rather  they  are  of  a 
mixed  wave  /  vortex  character  deriving  from  the  Tollmien-  Schlichting  eigenmode 
spectrum. 
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Linear  approximation 


The  theory  exposed  describes  essentially  nonlinear  motion.  When  turning  to 
stability  problems,  we  may  linearize  the  system  of  governing  equations  by  putting 

(u -  y, v,  w, p,d)  =  a («,  v,  (23) 


and  applying  the  Laplace-double  Fourier  transform 

[«c(  ® > hy,  m) . 4,(0 ,  k,  m )] 

0°  00  00  ^  (24) 

=  \d>c  \(fc\e~^UM+m^^u[i,x,y,!k),.r,At,x,Z) 

-00  -00  0 

Upon  solving  the  system  of  resulting  ordinary  differential  equations  the  dispersion 
relation  (Ryzhov,  2001;  Ryzhov  &  Bogdanova-Ryzhova,  2002) 


4*  1* 

4>(Cl)=  g(A,/^A£,,e3KZ>),  fi  =  Ad** 


(25a, b) 


to  connect  the  frequency  co  with  two  wavenumbers  comes  into  play.  Here 
cL4i{Q), 


(D(Q)  = 


dY 


(26) 


is  a  standard  function  depending  on  the  first  derivative  dAiipl )/ dY  and  an  improper 
integral 


Ppl)  =  jjii  Y)dY 


(27) 


n 


of  the  Airy  function  di(Y) .  The  arguments  of 


Q=P 


1^  +  m2 


5 

Id 


k2  +  (l- Nfc)'  rrt 


—  +  8  K  D\ 


(28) 


are  expressed  through  both  wavenumbers  and  do  not  contain  the  frequency. 

Spectral  side  band 

The  dispersion  relation  (25  a, b)  is  cast  in  the  form  based  on  the  scaling  intrinsic  to 
the  Tollmien-Schlichting  wave  spectrum.  However,  the  two  terms  in  curled  brackets  on 
the  right-hand  side  of  (28)  are  of  different  nature.  While  the  first  term  is  associated  with 
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the  wave  motion  and  gives  leading-order  contribution,  the  second  term  originates  from 
the  curvature  effects  and  brings  c(e 3 )  -correction  in  the  Tollmien-Schlichting  range  of 

spectrum.  Evidently,  a  side  band  of  frequencies  and  wavenumbers  should  emerge  from 
the  interaction  of  the  wave  and  vortex  eigenmodes  where  they  would  exert  an  equal 
impact  on  the  disturbance  pattern.  This  requirement  defines  the  spectral  interval 

1  2  _JL  2 


0  =  e7(k.C^7(d  g=  Re  28(kZ^7cd  g»  0.25co  g 

(29a) 

£  3  1 

k-  87(kZ^7&g=  Re  56(k Z^7&g«  0.12 kG 

(29b) 

_3  1  A  1 

m-  e  7 (k.C^-7  mG  —  Re56^/}?/?^,*  IS^rriQ 

(29c) 

of  the  spiral-type  G  rtler  vortices  for  typical  transition  conditions  in  an  incompressible 
boundary  layer  with  Re  ~  5  •  1 05 . 

The  auxiliary  variable  Q.  and,  hence,  0(fl)  are  invariants  and  the  affine 
transformation  (29a-c)  but  Q  and  D  change  to 
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with  both  terms  in  curled  brackets  on  the  right-hand  side  of  (30a)  becoming  of  equal 
order  in  magnitude.  Thus,  in  accord  with  the  above  conjecture  a  low- frequency  side  band 
comes  about  from  the  wave  /  vortex  interaction. 


Group  velocities 


The  centrifugal  force  effects  on  the  overall  system  of  pulsations  can  be  perceived 
from  a  comparison  of  the  first  dispersion  curves  in  the  plane  of  complex  frequencies. 
Figure  4  drawn  with  positive  values  of  k  is  typical  of  the  subsonic  boundary  layer  on  a 
flat  plate.  The  first  dispersion  curve  consists  of  a  single  lobe  starting  at  infinity  as  k-)  oo 
and  terminating  in  the  origin  with  k->  0.  The  shape  of  the  first  dispersion  curves  in  the 
plane  of  complex  G  rtler  frequencies  shown  in  Figure  5  for  Wq  =  0.5  and  three  values  of 

the  Mach  number  is  different.  The  formation  of  two  lobes  connected  through  a  small  loop 
on  top  is  a  common  feature  of  all  three  plots.  On  theoretical  grounds,  the  prediction  of  the 
behavior  of  disturbances  associated  with  the  unusual  two-lobe  shape  of  the  first 
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dispersion  curve  leans  upon  the  general  physical  concept  of  the  group  velocity  (Landau  & 
Lifshitz,  1959).  The  location  of  the  global  and  local  maxima  of  the  real  part  91  (co  of 
the  complex  frequency  in  the  co  G-plane  is  of  special  significance.  The  global  maximum 

of  91  (co  on  the  right-hand  lobe,  a  counterpart  of  that  in  Figure  4,  is  responsible  for  the 

generation  of  short-scaled  oscillation  cycles  in  the  wave  packet  sweeping  downstream. 
The  wave  packets  of  this  type  underlie  the  conventional  road  to  transition  provoked  by 

convective  instability.  A  tiny  positive  peak  of  91  (a>  q)  in  the  upper  portion  of  the  loops 

in  the  shape  of  the  first  dispersion  curves  in  Figure  5  is  endowed  with  a  negative  value  of 
the  group  velocity.  As  a  consequence,  this  peak  induces  a  wave  packet  consisting  of 
much  longer  oscillation  cycles  and  capable  of  advancing  upstream  of  a  site  where  the 
disturbance  was  given  birth. 

Absolute  instability  in  the  streamwise  direction 

Thus,  the  upstream  moving  modulated  signals  are  brought  about  by  the  interaction 
of  the  Tollmien-Schlichting  wave  and  G  rtler  vortex  eigenmodes.  Centrifugal  forces 
create  vortical  disturbance  structures  elongated,  in  keeping  with  (29b, c),  in  the 
streamwise  direction.  The  surface  curvature  provides  support  to  centrifugal  forces  which 
become  zero  if  a  flat  plate  is  under  consideration.  The  Blasius  boundary  layer  mostly 
studied  both  theoretically  and  experimentally  so  far  experiences  no  wave  /  vortex 
interaction  and  shows  no  absolute  instability  in  the  streamwise  direction.  However,  the 
first  evidence  of  upstream  influence  exerted  by  the  vortical  disturbances  from  a  concave 
insert  mounted  on  an  otherwise  flat  plate  can  be  traced  back  to  remarkable  wind-tunnel 
tests  by  Mangalam,  Dagenhart,  Hepner  &  Meyers  (1985). 

One  feature  of  the  first  dispersion  curves  in  Figure  5  is  worth  noting  beforehand 

to  obtain  reliable  results  from  computations.  The  size  of  tiny  positive  peaks  of  91  (w  q,) 
in  the  upper  portion  of  the  loops  in  the  shape  of  these  curves  is  much  less  compared  to  the 
magnitude  of  the  global  maximum  of  91  (co  positioned  on  the  right-hand  lobes.  Hence 

a  need  for  filtering  out  the  high-frequency  oscillation  cycles  becomes  obvious  if  both  the 
downstream  moving  highly-modulated  signal  and  the  upstream-advancing  wave  packet 
are  intended  to  be  computed  simultaneously.  The  plots  in  Figures  6,7  are  drawn  using  this 
numerical  technique.  They  correspond  to  two  dispersion  curves  in  Figure  5,  one  for  an 
incompressible  boundary  layer  and  the  other  for  a  boundary  layer  in  the  high  subsonic 
Mach  number  regime.  The  upstream  propagating  pulsations  predicted  above  theoretically 
solely  from  the  notion  of  the  negative  group  velocity  are  clearly  seen  in  both  plots.  The 
frequency  of  pulsations  grows  with  the  Mach  number  increasing.  The  absolute  instability 
in  the  streamwise  direction  may  provoke  earlier  transition  to  turbulence. 
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Conclusions 


The  present  study  leads  to  the  following  conclusions: 

•  Unsteady  spiral-type  G  rtler  vortices  make  up  a  low-frequency  side  band  in  the 
Tollmien-Schlichting  wave  spectrum. 

•  Absolute  instability  in  the  streamwise  direction  provoked  by  the  wave  7  vortex 
interaction  is  inherent  in  any  low  Mach  number  as  well  as  high  Mach  number 
subsonic  boundary  layers. 

•  Filtering  out  the  high-frequency  oscillation  cycles  is  required  to  compute 
simultaneously  the  wave  packets  moving  downstream  as  well  as  advancing  upstream, 
against  the  oncoming  flow. 

Future  work:  thermal  problem 

Temperature  disturbances  remained  mostly  beyond  the  scope  of  traditional 
hydrodynamic  stability.  However,  understanding  the  reasons  for  the  heat  transfer 
coefficient  to  increase  is  of  prime  importance  in  the  turbine  blade  design.  This  provides  a 
strong  impetus  to  extent  the  stability  analysis  to  thermal  wave  packets  propagating 
upstream.  The  complex-valued  amplitude  of  the  forced  temperature  pulsations  introduced 
in  much  the  same  way  as  in  (23)  and  (24)  obeys  an  inhomogeneous  Airy  equation 


j2  rr  t  2  2 

subject  to  the  limit  condition 
Tc  ->  ^  as  |lf  — >  o° 


(31) 


(32) 


at  the  upper  reaches  of  the  viscous  near-wall  sublayer  and  the  boundary  condition 
Tc  =  0  at  Y=  Q  (33) 


at  the  solid  surface  held  at  constant  temperature.  Here  vc  stands  for  the  normal-to-wall 
velocity  component  coming  from  the  above  solution  for  the  viscous  near-wall  sublayer.  A 
transformation  to  the  G  rtler  spectral  range 

-  i  _  I  _ 

Vc  =  6  7  (k47  VcG,  Tc  =  s  7  (k 4' 7  (34) 

preserves  the  boundary- value  problem  intact. 

Some  recent  results  based  on  a  solution  to  this  problem  are  submitted  by 
Bogdanova-Ryzhova  &  Ryzhov  to  ASME  Turbo  EXPO  2003.  As  it  follows  from  (31), 
the  temperature  oscillations  in  the  upstream  moving  wave  packets  may  have  even  larger 
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amplitude  than  the  normal-to-wall  velocity  (and  pressure)  if  the  length  of  G  rtler  vortices 
grows  without  bounds  as  kG~)  0. 
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Cooled  HP  turbine  rotor  blade 
(courtesy  of  Rolls-Royce  pic,  from 
S.L.  Dixon,  1998). 


Figure  1.  Typical  shape  of  a  turbine  blade. 


Sketch  of  the  streamwise  vortices  developing  on  a 
concave  wall  due  to  the  Gortler  instability  mechanism. 


Figure  2.  Steady  vortical  disturbances  in  the  classical  boundary  layer. 


Figure  3.  Unsteady  vortical  disturbances  in  an  interactive  boundary  layer. 


Figure  4.  Complex  frequencies  for  incompressible  boundary  layer. 
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m  =0.5 


3J(co) 

Figure  5.  Complex  frequencies  for  compressible  boundary  layers  in  different  subsonic 
regimes. 


m  =  0.5  M„=0.00 


Figure  6.  Downstream  and  upstream  propagating  disturbances  in  an  incompressible 
boundary  layer. 


m  =0.5  M -,=0.98 


t=13.00 


Figure  7.  Downstream  and  upstream  propagating  disturbances  in  compressible  boundary 
layer  in  the  high  subsonic  Mach  number  regime. 
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